Tailoring the Phonon Band Structure in Binary Colloidal Mixtures 
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We analyze the phonon spectra of periodic structures formed by two-dimensional mixtures of 
dipolar colloidal particles. These mixtures display an enormous variety of complex ordered config- 
urations [J. Fornleitner et a/., Soft Matter 4, 480 (2008)], allowing for the systematic investigation 
of the ensuing phonon spectra and the control of phononic gaps. We show how the shape of the 
phonon bands and the number and width of the phonon gaps can be controlled by changing the 
susceptibility ratio, the concentration and the mass ratio between the two components. 

PACS numbers: 63.20.Dj, 63.22.+m, 82.70.Dd 



Materials with a band gap in their spectrum of trans- 
mitted sound waves have been the focus of intensive re- 
search recently. Most often, these "phononic crystals", 
termed in analogy to the more familiar photonic crys- 
tals, are analyzed and fabricated on the basis of macro- 
scopic approaches: following methods from scattering 
theory, materials with a periodic modulation in their 
elastic properties and/or density are assembled, so that 
mismatches in the speed of sound and destructive inter- 
ference lead to the desired gaps [1-3]. An alternative 
and very promising route to understand and engineer 
phononic crystals is to focus on the microscopic prop- 
erties that directly shape the dispersion curves, i.e., the 
interparticle interactions. The latter uniquely determine 
both the ground-state configuration of the system and its 
elemantary excitations above the same, i.e., the phonons. 
A class of materials where this microscopic approach is 
particularly fruitful are colloidal crystals. In such sys- 
tems, the interactions are tunable and versatile. It has 
been shown that the dispersion curves of two-dimensional 
colloidal crystals can be shaped and controlled by suit- 
able external substrate potentials 0-[6|. Here, we fo- 
cus on systems that exhibit phononic gaps without the 
presence of external fields, namely mixtures: the non- 
trivial unit cells of their ground-states give rise to optical 
branches in the phonon band structure (PBS), thus open- 
ing the possibility to induce gaps by suitable adjustments 
in the interactions. Here, we explore the possibilities to 
tune the PBS's of a binary mixture of dipolar colloids via 
changes in the susceptibility ratio and concentration. 

Binary dipolar monolayers are readily available to ex- 
periments [7H12| and exhibit a rich wealth of stable or- 
dered structures [l3|, Q • Experimental realizations em- 
ploy super-paramagnetic particles of different suscepti- 
bilities Xi which are trapped at the interface of a pen- 
dant water droplet to ensure a planar geometry [12] . An 
external magnetic field B applied perpendicular to the 
water-air-interface induces magnetic moments in the col- 



loidal particles parallel to the external field, M$ = x^B, 
with i — A,B labeling the particle species. We empha- 
size that here the external field serves only as a means to 
influence the two-body interactions in the Hamiltionian 
and does not act on the systems as a one-body poten- 
tial, as is the case of laser beams in Refs. [4-6]. For the 
sake of simplicity, we assume that the physical size of 
the two particle species is the same and given by their 
common diameter a, the disparity in their interactions 
arising from different degrees of doping with ferromag- 
netic material, such as Fe203. Setting rrii = Xi/XA < 1, 
the ideal dipole-dipole repulsion acting in the mixture 
can be written as &ij(x) = smirrij/x 3 , i,j = A, £?, with 
x = r/cr and e = ^oXa!-^ 2 /^ 71 "^ 3 )- The ground state 
is thus determined by the asymmetry in dipole strength, 
given by m = Xb/xa < 1, and by the composition of 
the mixture, C = n^/^ + n#), with tia(b) being the 
number of strong (weak) dipoles per unit cell. 

Colloids in solution obey Langevin dynamics and their 
full equations of motion include the interparticle forces, 
the random collisions with the solvent and hydrodynamic 
interactions [TBI, involving the masses of the particles, 
the frictions constants and the solvent properties via the 
Navier- Stokes equations. Accordingly, the full dynamics 
of colloidal crystals have been thoroughly analyzed by 
theory and experiment in the past !l5l Il8j . Here we focus 
only on the excitation spectrum of the crystal, which is 
expressed in terms of the dispersion relation A(q), where 
q denotes a wavevector in the first Brillouin zone and 
A(q) the eigenvalue. Although oscillations in a colloidal 
crystal are overdamped, the quantities A(q) are still ex- 
perimentally measurable by means of the equipartition 
theorem [4| , a strategy that has been successfully applied 
to the one-component version of the system at hand [T(j • 

Accordingly, we apply harmonic lattice theory [Tj| 
|20| and obtain the dispersion curves A(q) by solving 
the eigenvalue equation A(q)c^ = Y^u'i' ^viyi'^Cv'i', 
where D v %yi'{c\) is the Fourier transform of the dynam- 
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ical matrix D V iyi'(n,n'), defined as: 



D uiyU 'i'(n,n f ) = 



d 2 $(r N ) 

du nv idu n ' V 'i' 



(i) 



Here, the index n is used to label the unit cell in the 
lattice and v runs over all particles within one unit cell, 
v = 1, riA+riB, u ni/ being the displacement of the v-Xh 
particle in the n-th unit cell of the lattice from its equilib- 
rium position. Finally, the index i denotes the Cartesian 
components, i = x,y, and $(r JV ) is the total potential en- 
ergy of the crystal at the harmonic approximation. The 
dimensionless eigenvalues X(q)(j 2 /e are determined for 
q- values along those axes of the first Brillouin zone that 
link the points of high symmetry within the same. 

In addition to affecting the ground states, the suscep- 
tibility ratio m can be used to fine-tune the appearance 
of a PBS corresponding to a given ground state, as the 
stable configurations are robust against small variations 
in m. In the following discussion, we focus on mixtures 
with C > 1/2, as the complexity in the emerging patterns 
increases with the number of weak dipoles in the system 
and the PBS's reflect recurrences in the ground- states. 

We start by demonstrating the influence of the suscep- 
tibility ratio m on the PBS for mixtures with a concen- 
tration of C = 1/2. We vary the susceptibility ratio in 
a range from m = 0.003 to m = 0.41 and calculate the 
PBS's for the ground states predicted in Refs. 13|, 
In Fig. [1] we show the phonon spectra for C = 1/2 
and increasing susceptibility ratio from top to bottom 
(m = 0.018, 0.04, 0.18, and 0.41). The corresponding 
ground states, together with their first Brillouin zones 
and the path along which the PBS's were determined, 
are depicted in Fig. [T] as well. At the chosen concentra- 
tion, the dipolar mixture develops three different types 
of ground states depending on m. For very low m- values, 
the stronger dipoles form a hexagonal lattice unaffected 
by the presence of the weaker species, the latter arrang- 
ing themselves in lanes and occupying interstitial sites. 
The associated PBS is characterized by a nonuniform dis- 
tribution of dispersion curves: the topmost four branches 
span a much broader eigenvalue range than the remain- 
ing four modes at the bottom of the band structure. 
In addition, the 'unfurled' upper branches are separated 
from the compressed lower region by a distinctive gap, 
see Fig. [TJa). Small increases in the susceptibility ratio 
m leave the ground state unchanged, despite of the fact 
that interactions involving the weaker species intensify 
[l3L [HI- I n the PBS, a small increase in m is reflected 
by a widening of the gap and a gradual expansion of 
the lower dispersion curves. Above a certain threshold 
in m, the weak dipoles are able to distort the hexagonal 
pattern of their stronger counterparts. This decrease in 
symmetry in the ground state is also reflected in the PBS, 
see Fig. QJb). At susceptibility ratios between m = 0.06 
and 0.18, a regular square lattice is predicted to be the 
energetically most favorable arrangement of dipoles [14 1. 




FIG. 1: Phonon band structures for a binary mixture of dipo- 
lar colloids with a concentration of C = 1/2 and susceptibility 
ratios of (a) m = 0.018, (b) m = 0.04, (c) m = 0.18, and (d) 
m — 0.28. Band gaps are marked by shaded regions. The 
corresponding ground-state is depicted on the top left of each 
band structure, with the species of stronger (weaker) dipoles 
given by yellow (red) spheres and two lattice vectors mark- 
ing the chosen unit cell. The sphere size reflects the dipole 
strength. On the bottom left of each band structure, a cell of 
the reciprocal lattice with its first Brillouin zone, marked by 
broken lines, is shown. The path of q- values along which the 
band structures were determined, is indicated by a full line. 
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In the PBS's associated to the quadratic configuration, 
the four dispersion curves cover the eigenvalue range in 
a uniform fashion: all branches are fully unfurled and 
no gaps occur. Variations in the susceptibility ratio that 
leave the ground-state unchanged result in a stretching 
of the PBS without altering its general appearance. As 
an example, Fig.[TJc) shows the PBS of the square lattice 
formed at m = 0.18. For m > 0.28, the ground-state of 
the binary mixture changes from the regular square lat- 



tice to the ^-structure [2l|, which consists of distorted 



hexagons of strong dipoles accommodating two particles 
of the weaker species. The PBS's of the ^-lattice are 
again characterized by the occurrence of distinct eigen- 
value bands. The number of these bands and the width 
of the gaps separating them from each other are tunable 
via the susceptibility ratio: at m = 0.28, i.e., close to the 
transition from the square to the ^-lattice, three gaps 
are observed, see Fig. [if d). 

At higher concentrations of weak dipoles, more exotic 
particle arrangements with a low degree of symmetry ap- 
pear among the ground states. For the PBS's, the trends 
described above become more pronounced. In Fig. [2] we 
show selected PBS's obtained for mixtures with (7 = 7/9 
and susceptibility ratios in the range from m = 0.003 to 
m = 0.28. At small values of m, the non-uniform distri- 
bution of dispersion curves is recovered. The four curves 
at the top of the PBS again cover a much broader eigen- 
value range than the rest of the branches and a gap sep- 
arating the compact from the extended dispersion curves 
is observed, even at very low values of m. The suscep- 
tibility ratio can again be used to tune the gap width, 
see Figs. [D(a), [2fb). In contrast to the previous case, 
C = 1/2, the distortion of the underlying hexagonal lat- 
tice induced by an increase of m, does not destroy the 
separation of the PBS into distinct bands. A further in- 
crease in m opens, on the contrary, additional gaps, first, 
by causing branches from the lower, compressed region to 
detach and then, as m increases, by a progressive flatten- 
ing of the topmost curves, going from Figs. [2fc) to[2fd). 
At m = 0.09, Fig.[2fd), a maximum of five distinct bands 
is reached. If the susceptibility ratio is increased above 
a certain threshold, i.e. m > 0.18, the ground state of 
the mixture undergoes a transition to lane-like arrange- 
ments, forcing all gaps in the PBS to close and leading 
to a uniform distribution of well modulated dispersion 
curves, Fig. [2fe). 

The general trends become even clearer in the sequence 
of PBS's obtained at C = 4/5, see Fig. H for PBS cal- 
culated at m = 0.003, 0.018, 0.09, and 0.18. The band 
structures at low values of m are characterized by two 
distinct bands, the width of the gap separating them is 
well tunable via the susceptibility ratio, Figs. [3fa), [3fb). 
At intermediate values of m, the dispersion curves flatten 
out and multiple gaps appear between the less modulated 
branches, see Figs. EUc), [3](d), which vanish again with 
the onset of lane formation observed at high susceptibil- 




FIG. 2: Phonon band structures for a binary mixture of dipo- 
lar colloids with a concentration of C = 7/9 and susceptibility 
ratios of (a) m = 0.003, (b) m = 0.018, (c) m = 0.04, (d) 
m 0.09, and (e) m 0.28. 
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FIG. 3: Phonon band structures for a binary mixture of dipo- 
lar colloids with a concentration of C = 4/5 and m- values: (a) 
m = 0.003, (b) m 0.018, (c) m 0.09, (d) m 0.18. 

ity ratios. 

Thus far, we have not considered a possible mass asym- 
metry of the two species in our calculations. This is done 
because the colloidal systems that allow for an experi- 
mental realization of our findings are Brownian systems, 
in which the mass of the particles is rendered irrelevant 
by overdamping 0, [lcj. Including the particle masses 
M v by multiplying expression (pQ) for the dynamical ma- 
trix with 1 / y/M v M v >, we can easily switch to Newtonian 
systems in our theoretical approach and gain the mass 
ratio M = Mb /Ma as an additional parameter to tune 



the PBS's with. In ongoing investigations we have found 
that mass asymmetry (as it occurs, for instance, in dusty 
plasmas) enhances the effects described above: by lower- 
ing M, the topmost 2ua branches can be decoupled from 
the rest of the band structure and shifted to higher eigen- 
values, thus opening a gap of variable width. In addition, 
the flattening of dispersion curves gets more pronounced 
with increasing mass- asymmetry, so that even completely 
flat branches can be induced on the PBS's. 

In summary, we have investigated the possibility to 
control the PBS's of a dipolar binary mixture via the pa- 
rameters determining the interactions and the ground- 
state of the system. By calculation of the dispersion 
curves for the stable particle arrangements, we could 
show that, for systems with a majority of weaker dipoles, 
multiple gaps can be opened by suitable adjustments in 
the susceptibility ratio m. In addition, remarkably flat 
branches were observed in the PBS featuring multiple 
gaps. Our results are easily verifiable in experimental se- 
tups that employ two-dimensional binary magnetic col- 
loids, where dispersion curves can be measured by means 
of the equipartition theorem 0, [To| . 
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